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Abstract 

We study the mass dependence of various quantities (like the average and maximum 
density, colhsion rate, participant-spectator matter, temperature as well as time zones for 
higher density) by simulating the reactions at the energy of vanishing flow. This study is 
carried out within the framework of Quantum Molecular Dynamics model. Our findings 
clearly indicate an existence of a power law in all the above quantities calculated at the 
balance energy. The only significant mass dependence was obtained for the temperature 
reached in the central sphere. All other quantities are rather either insensitive or depend 
weakly on the system size at balance energy. The time zone for higher density as well as 
the time of maximal density and collision rate follow a power law inverse to the energy of 
vanishing flow. 



1 Introduction 

It is now well established that the interactions at low incident energies are dominated by the 
attractive part of the nuclear mean field causing the emission of particles in the backward 
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angles. These interactions, however, become repulsive at higher incident energies that pushes 
the particles in the forward (positive) angles. While going from the low incident energies to 
higher incident energies, there is a particular energy at which the net flow disappears PP. This 
energy, (termed as "balance energy") Ebai has been found to be of significant importance for 
the understanding of the nature of nuclear interactions and related dynamics [2-24]. 

Recently, the Ehai was measured in ^^'^Au+^^''Au collisions [HIH], extending the mass range 
of Ehai between 24 and 394 units. In addition to ^^^Au+^^''Au coUisions, one has also measured 
the Ebai in '^C+^^C fl], 20Ne+27Al 7], ^^Ar+^^Al |1[I2], ^"Ar+^^AIE], ^^Ar+^^Sc 0[7l[I3], 
40Ar+5iY 642n+27Ai ^n^^ ^OA^+ssNi ||, 642^+48^^ ^ 58Ni+58Ni 001131, ^^Fe+^^Fe 
Hg, '^^Zn+ssNi ^'^Kr+^^Nb 0171, ^^Nb+^^Nb 0, ^^^Xe+'^^Sn and ^^^La+^^^La ^ 
collisions. Most of the above reactions are symmetric and central in nature. Some attempts are 
also reported in the literature that deal with the impact parameter dependence of the balance 
energy E^ai |SlinilHllIIllI3IIZll2niini- 

Interestingly, most of the theoretical attempts for disappearance of flow are within the 
Boltzmann-Uehhng-Uhlenbeck (BUU) model PEHilSllZllHllIIlliaiiaiinillllllHl- Some 
attempts, however, are also made within the framework of Quantum Molecular Dynamics 
(QMD) model |I3l IHl 1201 IHl 122 123 121 ■ Note that among all these attempts, only a few 
deal with the mass dependence of the disappearance of flow IH EllZl (HI 1221 123 1211 ■ There one 
reported a power law behavior (ocA"^) in the Efyai- For the first time, a complete study of the 
mass dependence of balance energy was presented by us where as many as sixteen systems, 
with mass between 47 and 476 were analyzed [221 123 121 ■ The excellent agreement between 
the experimental measurements and theoretical calculations allowed us to predict the balance 
energy in 238u_j_238u collision around 37-39 MeV/nucleon [221 • None of the above studies was 
extended to other heavy-ion phenomena at balance energy. Refs. [TH IT^ ITHl (23 1^ , give some 
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information about the nature of other variables at the balance energy. 

We here present a complete analysis of the nuclear dynamics at balance energy by analyzing 
more than 14 (nearly symmetric) reactions as reported in ref. [221 121] • Our present motivation is 
to investigate whether other dynamical variables (apart from the disappearance of flow) show a 
mass dependence at the balance energy or not. This present study is made within the framework 
of Quantum Molecular Dynamics (QMD) model [inil2nil2Ill22ll2nil21ll2nil2nil23l2Hll2i- 
Section II describes the model in brief. Our results along with the discussion are presented in 
section III. We summarize the results in section IV. 



2 The Model 

The present study is made within the framework of the Quantum Molecular Dynamics (QMD) 
model [211 123 121] where each nucleon interacts via two- and three-body interactions that pre- 
serve the nucleon-nucleon (nn) correlations and fluctuations. Here each nucleon is represented 
by a Gaussian wave packet with width ^J~L centered around the mean position rj(t) and the 
mean momentum pj(t): 

0.(r,P,t) = ^_l^e[-(^^--W)^He[-W^ (1) 
The Wigner distribution of a system with {^At^Ap) nucleons is given by: 

/(P.p.*) = 'ff'^el-'-'OIV-lx 



g[-lp-ftWr^Wft (2) 



The nucleons propagate under the classical equations of motion: 

dvi _ d{H) _ 
dt dpi ' 



(3) 



(ipi_ ^ d{H) 

dt dr. ' ^ ^ 



The Hamiltonian (H) is given by: 



{H) = {T) + {V); 



with 

^total ^local _j_ ^Yuk _|_ -^Coul /g\ 

ij ij '^J ^ ' 

Here 1^^°^^"', V^^°"' and Vj^^'^ stand, respectively, for the Skyrme, Coulomb and Yukawa parts 
of the nn interaction. Following refs. [221 1231 1211 ' ^ hard equation of state along with 

energy independent cross-section of 40 mb strength. This combination is reported to explain 
the experimental balance energy nicely \2'6\ . 



3 Results and Discussion 

As stated in ref . [231 , a hard equation of state along with energy independent nn cross-section 
of 40 mb strength yields a power law behavior ocA^. The power law {c.A^) over the exper- 
imental points yields r = —0.42079 ± 0.04594, whereas our theoretical calculations with nn 
cross-section of 40 mb strength had r = —0.41540 ± 0.08166 (231 • It is worth mentioning that 
this was the closest agreement obtained so far. For the present mass dependent analysis, we sim- 
ulated the reactions of ^^Ne+^^Al (b=2.6103 fm), ^^Ar+^^Al (b=2 fm), ^"Ar+^^Al (b=1.6 fm), 
40Ar+45gp (b=3.187fm), ^^Ar+^iV (b=2.442 fm), ^OAr+^^Ni (b=0-3 fm), ^^Zn+^^rji (^^2 fm), 
58Ni+58Nj (b=2.48 fm), ^^Zn+^^m (b=2 fm), s^Kr+^^Nb (b=4.07 fm), ^^Nb+^^^b (b=3.104 
fm), i29Xe+^'^*Sn (b=0-3 fm), ^^^La+^^^La (b=3.549 fm) and i97Au+i97Au (b=2.5 fm) at their 
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corresponding theoretical balance energy^ which is, respectively, 119, 74, 67.3, 89.4, 67.8, 64.6, 
59.3, 62.6, 56.6, 69.2, 57, 49, 51.6, and 43 MeV/nucleon. The reactions were followed till nuclear 
transverse flow saturates. As noted from the above, the balance energy is smaller in heavier 
colliding nuclei, compared to lighter one. As a result, one would expect early saturation in 
lighter colliding nuclei compared to heavy one. Though, the transverse flow saturates much 
early in lighter nuclei, some of the quantities, however, keep changing, therefore, we follow all 
the reactions uniformly till 300 fm/c. 

In the following, we shall first study the time evolution and then shall present the mass 
dependence of different quantities. 

3.1 The Time Evolution 

One of the motivations behind studying a heavy-ion collision is to extract the information 
regarding the hot and dense nuclear matter. In our approach, the matter density is calculated 
by H; 

/'('',*) = 'E"p^e<-«'-<'))V-). (7) 

Here At and Ap stand, respectively, for the target and projectile. In actual calculations, we take 

a sphere of 2 fm radius around the center-of-mass and compute the density at each time step 

during the reaction using eq. (7). Naturally, one can either extract an average density (p""^) 

over the whole sphere or a maximal value of the density {^p"*"^) reached anywhere in the sphere. 

In fig. 1(a), we display the (p"^^)/po whereas fig. 1(b) shows the (p™'"^)/po as a function of the 

reaction time. The displayed reactions are ^ONe+^^Al (A=47), ^^Ar+^^Sc (A=85), ^^Zn+^^Ni 

(A=122), 93Nb+93Nb (A=186), ^^^La+^^^La (A=278) and ^^^Au+^^^Au (A=394) spreading 
^The theoretical balance energy was calculated by extrapolating the flow at two different energies with a 

step of ±10 MeV/nucleon 



5 



over the whole mass range. We see that the maximal p""^^ is slightly higher for lighter systems 
compared to heavy ones. A similar trend can also be seen for the evolution of p^°-^ . Further, the 
maximal value of the density for medium and heavy systems is comparable with the average one. 
This clearly indicates that the dense matter is formed widely and uniformly in the central region 
of 2 fm radius. On the other hand, substantial difference in two densities for lighter colliding 
nuclei indicate that the dense matter is not homogenous and uniform in these reactions. Due to 
high incident energy, ^''Ne+^^Al reaction finishes much early compared to ^^^Au+^^''Au which 
is simulated at a relatively lower incident energy. Similarly, the peaks in (the maximum l^p^°-^^ 
and average (p"^^)) densities are also delayed in heavier nuclei compared to lighter one. The 
spreading of the high density zone in heavier colliding nuclei over the long time span indicates 
the on going interactions among nucleons. This result is in agreement with |29j . 

Another quantity directly linked with the density is the collision rate. In fig. 2, we display 
the dNcoii/dt versus time. This rate represents the net collisions after fulfilling the Pauli prin- 
ciple. Naturally, the attempted rate will be much higher than the allowed one. Due to larger 
interaction volume, the interactions among nucleons in heavy nuclei continue for a long time. 
This effect should be obvious if one looks the density profile (see fig. 1). A finite density zone 
will naturally lead to more and more nn collisions and as a result, the collision rate will be 
more for heavy colliding nuclei. 

As stated above, all the reactions are simulated at the balance energy where attractive forces 
counter balance the repulsive forces. This fact should also be reflected in the quantities like 
the spectator and participant matter. All nucleons having experienced at least one collision are 
counted as participant matter. The remaining matter is labeled as the spectator matter. The 
nucleons with more than one collision are labeled as super-participant matter. This definition 
gives us possibility to analyze the reaction in terms of participant-spectator fireball model. 
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In fig. 3, we display the normalized spectator matter (upper part) and participant matter 
(lower part) as a function of the time. At the start of the reaction, all nucleons constitute 
spectator matter. Therefore, no participant matter exists at t=0 fm/c. Since the ^'^Ne+^'^'Al 
reaction happens at a relative higher energy (=119 MeV/nucleon), the transition from the spec- 
tator to participant matter is swift and fast. On the other hand, due to low bombarding energy, 
the transition from the spectator to participant matter in heavier colliding nuclei is rather slow 
and gradual. Interestingly, at the end, all reactions (that happen between incident energy 43 
and 119 MeV/nucleon) lead to nearly same participant matter indicating the universality in 
the balancing of attractive and repulsive forces. 

From the above facts, it is clear that heavier colliding nuclei (at Ehai) have delayed and ex- 
panded evolution of the density and participant-spectator matter. It will be of more interest to 
see how their mass dependence behaves like. This will be discussed in the following paragraphs. 



3.2 The Mass Dependence 

In fig. 4, we display the maximal value of {p""^^) and l^p^"-^^ versus composite mass of the 



At~Ai 



< 



system. Note that all the reactions, considered here, are symmetric i.e. ri 

0.2. Very interesting, we see that the maximal value follows a power law (xA'^ with r being 

—0.05182 ± 0.00776 for average density (p"^^) and —0.11477 ± 0.01217 for maximum density 

^pmaxy2 jj-^ other words, a slight decrease in the density occurs with the size of the system. 

This decrease is much smaller compared to the Et,ai ijexpt = —0.42079 ± 0.04594 and Tth = 

—0.41540 ± 0.08166). Had these reactions being simulated at a fixed incident energy, the trend 

would have been totally different [21]. Since lighter nuclei cannot be compressed easily, their 

maximal density at a fixed incident energy will be less compared to heavy nuclei. Since Ebai in 
small deviation can be seen in the cases where f] 0. 
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the present case scales as A"^, a weak mass dependence is also observed in density profiles. 

The mass dependence of the (allowed) nn collisions is depicted in fig. 5. Here one sees a 
(nearly) linear enhancement in the nn collisions with the size of the interacting system. This 
enhancement can be explained mainly using a power law ocA^; with r = 0.87829 ± 0.01833. If 
one keeps the incident energy fixed (e.g. in the figure, we kept E=50 MeV/nucleon), the nn 
collisions should scale as "A". Our fitting gives r = 1.04633 ± 0.01712, which is very close to 
unity. 

A dynamical quantity that can serve as an indicator of the role of repulsive and attractive 
forces is the participant and spectator matter. Naturally, the possibility of a collision will 
depend upon the mean free path of nucleons. Similar is the case of spectator and participant 
matter. In fig. 6, we display the spectator, participant and super-participant matter (obtained 
at 300 fm/c) as a function of the total mass of the system. Very interesting, we see a nearly mass 
independent behavior of the participant matter (r = —0.03621 ± 0.00954). Similar behavior 
occurs in the case of spectator matter (r = 0.08323 ± 0.02232). A slight deviation occurs in 
the case of ^°Ne-|-^''Al. Some small fluctuations can also be due to the variation in the impact 
parameter, which is not fixed in the present study. The choice of impact parameter is guided by 
the experimental measurements. As noted in ref. [21], the variation in the impact parameter 
can have drastic influence on the participant/spectator matter. The super-participant matter 
shows a little more dependence on the mass of the system (r = —0.14241 ± 0.04184). This can 
be understood again by looking at the density profile (fig. 4). There one concluded that the 
lighter nuclei lead to higher densities. In other words, the mean free path will be smaller in 
lighter nuclei that results in more nn collisions. One should keep in the mind that the mass 
independent nature of the participant matter is not a trivial observation. For a fixed system 
mass, the participant matter depends linearly on the incident energy. In the present case, 
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though the mass of the system increases, their corresponding energy decreases, resulting in the 
net mass independent nature. This also suggests that the repulsive and attractive forces at 
Ehai counter-balance each other in such a manner that the net participant matter remains the 
same. One may also say that since the contribution of the mean field towards transverse flow 
is nearly mass independent [231201, one needs same amount of participant matter to counter 
balance the attractive forces. In other words, the participant matter can act as a barometer to 
study the balance energy in heavy-ion collisions. 

The associated quantity linked with the dense matter is the temperature. In principle, a 
true temperature can be defined only for a thermalized and equilibrated matter. Since in a 
heavy-ion collision, the matter is non-equilibrated, one can not talk of "temperature". One 
can, however, look in terms of a local environment only. In our present case, we follow the 
description of the temperature given in refs. [211 1^ E2] ■ Several different authors have given 
different descriptions of the local or global temperatures [^21 OH I2H1 ESI- Some studies of 
temperature are based on the fireball model jSSj, whereas others take the hydrodynamical 
model into account [31]. In refs. |34[ I35j. the thermalization is directly connected with the 
non-diagonal elements of the stress tensor. One has even defined the "transverse" temperature 
in terms of (Pj/2m); Pf is the average transverse momentum squared [21]. In the present 
case, the extraction of the temperature "T" is based on the local density approximation i.e. 
one deduces the temperature in a volume element surrounding the position of each particle at 
a given time step [211 Ell I22|- Here, we postulate that each local volume element of nuclear 
matter in coordinate space and time has some "temperature" defined by the diffused edge of 
the deformed Fermi distribution consisting of two colliding Fermi spheres which is typical for 
a non-equilibrium momentum distribution in heavy-ion collisions. 

In this formalism, (dubbed as hot Thomas- Fermi approach [21]), one determines the exten- 

9 



sive quantities like the density, kinetic energy, as well as the entropy with the help of momentum 
distributions at a given temperature. For more details, reader is referred to refs. [2^1 IHH 1^ - 
Using this formalism, we also extracted the average and maximum temperature within a central 
sphere of 2 fm radius as described in the case of density. 

In fig. 7, we plot the maximal value of (T"''^) and l^T^"-^) as a function of the composite 
mass of the system. Some fluctuations can be seen in these plots that can be due to the 
choice of impact parameter as well as incident energy [211 E21 IHH IHZl IHH] • As stated above, 
the impact parameter choice is guided by the experimental constraints. Further, the E},ai was 
extracted using a straight line interpolation, therefore, both these factors may add to the 
present fluctuations. One sees that both quantities can be parameterized in terms of a power 
law function (xA^; The power factor r is rather quite large (being equal to —0.83743 ± 0.11355 
and —0.51079 ± 0.08218), respectively, for the average and maximal temperature. This sharp 
mass dependence in the temperature is rather in contradiction to the mild mass dependence 
obtained in all other quantities. This is not astonishing since temperature depends, crucially, 
on the kinetic energy (or the excitation energy) of the system (23 EH E21 EH EHl • It was shown 
in ref. j27| EH] that for a given colliding geometry, the maximal value of the temperature does 
not depend upon the size of the interacting source. Rather it depends only on the bombarding 
energy. 

In fig. 8, we display the time of maximal collision rate and average density. We see a 
power law behavior in both the quantities. The small balance energy in heavy nuclei delays 
the maximal compression. Interestingly, the power factor r is close to (1/3) in both the cases. 
The Ebai was shown to scale with power factor r —0.4. In other words, the time of maximal 
collision rate and density varies approximately as a inverse of the Ei,ai. 

Apart from the maximal quantities, another interesting quantity is the dense zone at the 
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balance energy. This is depicted in fig. 9 where we display the time interval for which p>po 
(upper part) and p>po/2. Again both quantities follow a power law behaviour. Interestingly, 
the time intervals for the high density have a power law dependence with r = 0.33853 ±0.05358 
and T = 0.46833 ± 0.04265, respectively, for p>po and p>po/2 which is again very close to the 
inverse of the mass dependence of -Em- This also points toward the fact that the formation 
and identification of the fragments is delayed in heavier nuclei compared to lighter nuclei. This 
conclusion is in agreement with earlier calculations I38j. 

4 Summary 

Using the QMD model, we presented the mass dependence of various quantities (such as aver- 
age and maximum central density, temperature, collision dynamics, participant and spectator 
matter as well as the time zone for hot and dense nuclear matter) at the energy of vanishing flow 
(Ehai). This study was conducted using a hard equation of state along with nn cross-section of 
40 mb strength. This combination is reported to explain the experimentally extracted balance 
energy for large number of cases [221 • calculations present several interesting facts: 

The reaction saturation time is smaller for lighter nuclei compared to heavy nuclei. The 
maximal values of density, temperature and collision rate is also shifted accordingly. In all the 
cases (i.e. in average and maximum central density, temperature, participant and spectator 
matter etc.), a power law dependence can be seen. The only quantity where power factor r is 
significant (with t>|0.2|) is the temperature reached in the central zone. Other quantities are 
nearly mass independent. The mass independent nature of the participant matter makes it a 
good alternate indicator for determining the balance energy. The existence of dense zone scales 
as inverse of the energy of vanishing flow. 
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Figure Captions: 



Fig. 1: The evolution of (a) average density {p"-""^) and (b) the maximum density l^p^"'^) 
reached in a central sphere of radius 2 fm as a function of time. Here reactions of ^'^Ne+^^Al 
(b=2.6103 fm), ^^Ar+^^g^ (b=3.187 fm), 64Zn+5«Ni (b=2 fm), ^SNb+^^^Nb (b=3.104 fm), 
^^^La+^^^La (b=3.549 fm) and ^^'''Au+^^^Au (b=2.5 fm) are simulated at their correspond- 
ing theoretical balance energies (for details, see the text). The shaded area represents the 
reaction of 2°Ne+^'''Al. 

Fig. 2: Same as fig. 1(a), but the rate of allowed coUisions dNcou/dt versus reaction time. 

Fig. 3: Same as fig 1, but the time evolution of normalized spectator matter (upper part) 
and participant matter (lower part). 

Fig. 4: The maximal value of the average density {p°''"^)max (upper part) and maximum 
density {p"^"'^)max (lower part) as a function of the composite mass of the system. The solid 
lines are the fits to the calculated results using c.A'^ obtained with minimization. 

Fig. 5: The total number of the allowed collisions (obtained at the final stage) versus 
composite mass of the system. The solid squares and open diamonds are the results obtained 
at Ehai and 50 MeV/nucleon, respectively. The sohd and dashed lines are the fits obtained with 
procedure explained in fig. 4. 

Fig. 6: Same as fig. 5, but for the final saturated participant, spectator and super- 
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participant matter. 

Fig. 7: Same as fig. 5, but for the maximal value of the average temperature (upper part) 
and maximum temperature (lower part). 

Fig. 8: Same as fig. 5, but for the time of maximal value of collision rate (open stars) and 
average density (solid squares). The solid and dashed hnes represent the fits with power law 

Fig. 9: Same as fig. 8, but for the time zone for {p>Pq) (upper part) and for {p>po/'2) 
(lower part) as a function of composite mass of the system. 
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Fig 1 : "Nuclear Dynamics at the Balance Energy" by Sood and Purl 
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Fig 2: "Nuclear Dynamics at the Balance Energy" by Sood and Puri 
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Fig 4: "Nuclear Dynamics at the Balance Energy" by Sood and Puri 
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Fig 5: "Nuclear Dynamics at the Balance Energy" by Sood and Purl 
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Fig 6: "Nuclear Dynamics at the Balance Energy" by Sood and Purl 



Fig 7: "Nuclear Dynamics at the Balance Energy" by Sood and Purl 
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Fig 9: "Nuclear Dynamics at the Balance Energy" by Sood and Purl 
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